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Abstract 

The classical approach to solvability of a mathematical problem is to define a method which 

includes certain rules of operation or algorithms. Then using the defined method, one can show 
that some problems are solvable or not solvable or undecidable depending on the particular 
method. With numerical solutions implemented in a computer, it might be more practical 
to define solvability of a mathematical problem as a complex decomposition problem. The 
decomposition breaks the data into a simple part and a complex part. The simple part is the 
solvable part by the method prescribed in the problem definition. The complex part is the 
leftover of the simple part. Complex part can be viewed as the "residual" of data or operator. 
It should be interpreted and not to be discarded as useless. We will give different examples 
to illustrate the more practical definition of solvability. The complex part is not noise and 
should not be viewed as useless part of the data. It has its own merit in terms of topological or 
geological interpretation. We have de-emphasized absolute solvability and have emphasized the 
practical solvability where the simple and complex parts both play important roles. 

Introduction 

We will start with some classical examples of solvability in mathematics. 

The simplest example is the unsolvability of square root of 2 with rational numbers. The 
prescribed method is to form fractions m/n where m and n arc integers and n is nonzero. Square 
root 2 is not solvable with this method of rational numbers or fractions. This is normally proved 
by way of contradiction. We can assume that square root 2 is a rational number. Then find a 
contradiction in the end. This is not a constructive proof. 

The next example is solving polynomials. From the Fundamental Theorem of Algebra, it is 
known that a polynomial with complex coefficients could be factored into linear parts. However, 
the theorem is an existence theorem which is not constructive but only shows the existence of the 
linear factors. 

Before computers, we had to solve the equation with only simple combinations. We restricted 
the definition of solvability of polynomials as solvable by radicals. This means that we can only use 
radicals (square root, cube roots, squares, cubes etc) to form the solution. For quadratic equation. 



we have the quardratic formula which solves the quadratic equation with square roots and powers 
of the coefficients. The quadratic is solvable by radicals. But the higher order polynomials are not 
solvable by radicals in general. This leads to group theory, more specifically Galois group theory 
of permutations, which gives a more abstract formalism of solvable by radicals. 

We then extend Galois group to more general group theory to explain physical problems as a 
group of transformations on physical equations. The group theoretical application can be viewed 
as a group of symmetries or a group of operators like orthogonal matrices. Then physics or other 
physical phenomenon is studied through the eyes of transformation groups. It has been a successful 
approach to many physical problems. 

The challenge to the group theoretic approach is that systems have to be described as reversible 
since each element of a group has an inverse. We know that not all matrices are invertible. For 
physical problems which are "lossy" , we have to introduce a more general formalism like semigroup. 
The semigroup of transformations allows noninvertible elements or operators. The collection of 
square matrices (e.g. all 2 x 2 matrices) form a semigroup under matrix multiplication but is not 
a group. Semigroup can be used to view the world as non-reversible processes. 

Unsolvability of Cantor layers in seismic imaging 

Before we look into different ways to express simplicity/ complexity, let us examine the Cantor layer 
construction in seismic imaging. It is taking the Cantor middle third construction of layering. See 
Lau and Yin SEG 2005 for more details. 

Figure 1 shows different stages of Cantor layers and the corresponding primary and multiple 
reflections. The Cantor layers converge to the Cantor set of reflections C. The multiple reflections 
converge to a continuum, i.e., 

C + C = I, 

where C is the Cantor set and / is a continuum. 

Cantor set is both sparse and dense. If the Cantor set is the support of the primary reflection 
coefficients, then C + C is the support of the multiple scattered refiections which is a continuum. 
Hence it is not possible to solve for Cantor coefficients which is sparse. The multiples form a 
continuum which is the opposite of sparsity. 

This is an example in which there is no solution. So practical approach would be needed to 
simplify the problem. 

Reduce data complexity (with predetermined basis) 

The simplest way to reduce data is to do a linear filter based on a predetermined basis like Fourier 
sinusoids or wavelet transform. The complex decomposition is to apply the bandpass filter to output 
the simple part. The difference of the input and output is the complex part. The predetermined 
basis has the advantage of being fast computationally. It has the disadvantage of imposing the 
geometry of the basis to the data. So if Fourier transform is used, the simple part would look like 
sinusoids. Geologically it is not likely to have sinusoids which imply recurrent geologic events over 
many millions years. 

Figure 2 shows 3 panels where the first panel is the input seismic image, the second is the output 
of the simple part of the seismic and the last panel is the difference which is the complex part of 
the data. 
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Figure 2: The first panel is the input seismic image, the second is the output of the simple part of 
the seismic and the last panel is the difference which is the complex part of the data. 
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Figure 3: The first panel is the input seismic image, the second is the output of the simple part of 
the seismic and the last panel is the difference which represents the complex part of the data. 

Reduce data complexity (without predetermined basis) 

The smooth part of the data can be thought of as the simple part of the data. We have used 
variational method to estimate the simple part of the data which can be used for further processing 
of seismic data. We have employed the formulation of Chan et al for minimizing the LI of the first 
derivative. 

Another approach is to use diffusion semigroup, e.g., Lau, Yin, Coifman and Vasilliou SEG 
2009. We will show just the variational method in this paper. 

Figure 3 shows 3 panels where the first panel is the input seismic image, the second is the output 
of the simple part of the seismic and the last panel is the difference which represents the complex 
part of the data. 

Even though the third panel of the complex part looks noisy, it is interpretable in terms of 
the details of the diffractions. The geometry of the diffractions can be interpreted and it could be 
traced through the diffraction limps. The actual amplitudes of the diffractions cannot be explained 
completely by a wave equation. This shows the geometric complexity of the diffractions. 

The seismic data is not solvable in an absolute sense. Not every sample can be explained by 
the wave equation. The simple part is what is solvable with the current algorithm. The complex 
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Figure 4: Seismic responses to migrate the recorded data to their proper position. The first panel 
is the correct operator generating the impulse response. The second panel is the simple operator 
using a single velocity using fk migration. The third panel is the complex part of the operator 
which is the difference of the first and second panel. 

part should not be discarded but should be interpreted. 

Reduce operator complexity (with predetermined basis) 

Let us start by introducing a simple semigroup which is just the exponential. We know that 
exp(t + s) = exp{t)exp{s). We can reduce wave equation from a second order partial differential 
equation into an exponential by performing Fourier transform on space and time domain. Then the 
depth information is obtained by an exponential phase shift in Fourier domain and then transform to 
depth domain. This is correct if it is simple velocity field like velocity field with no lateral variation. 
This is commonly called fk migration in geophysics which is using Fourier (/ for frequency and k 
for wavenumber) to solve the pde with predetermined basis like Fourier basis. Vertical variation is 
solved by recursive application of fk migration. 

Figure 4 shows 3 panels of seismic responses to migrate the recorded data to their proper 
position. The first panel is the correct operator generating the impulse response. The second panel 
is the simple operator using a single velocity using fk migration. The third panel is the complex 
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part of the operator which is the "difference" of the first and second panel. 

Reduce operator complexity (without predetermined basis) 

If we do not use a predetermined basis hke Fourier transform, we can use the velocity model 
coupled with the wave equation to form a simple operator. See Sandbcrg and Beylkin 2009. A less 
ambitious approach is to lump two variables into one, e.g., velocity and thickness can be lumped 
into single variable of travel time. See Lau and Yin SEG 2003 on " micro-tomography" , Lau and 
Yin arXiv 2010 on "Transformation Semigroup and Complex Topology". Lumping variables is a 
way to reduce operator complexity without using a predetermined basis. 

Reduce topological complexity 

One final approach is to reduce topological complexity. We use computational topology to define 
Bctti numbers. The different Betti numbers of different dimensions indicate the number of "holes" 
for that dimension. Sec Lau and Yin SEG 2009 on "topological simplicity" and arXiv 2010 on 
"L0-I-L1+L2 optimization". Minimizing the Betti numbers is a way to reduce topological com- 
plexity. By the same token, topologically complex data cannot be solved by wave equation in an 
absolute sense. 

Conclusion 

We have used a more practical definition of solvability. Absolute solvability is in general difficult to 
decide. A data set or an operator could be reduced in complexity by various means as demonstrated 

above. Complex decomposition breaks up the dataset or an operator into simple part (solvable part) 
and complex part (interpretable part). The complex part should be preserved for interpretation 
and should not be discarded as "noise". They have topological and geological meaning in terms of 
interpretation. 
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